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Summary: Some topological properties of stochastic flow (pt{x) generated by stochastic dif- 
ferential equation in a Ml with normal reflection at the boundary are investigated. Sobolev differ- 
00 . 

. entiability in initial condition is received. The absolute continuity of the measure-valued process 

^ ■ ^ o ip~[^ , where /i <C A"^, is studied. 

tj ' Flows generated by SDEs in Euclidean space is a well-studied topic nowa- 

O ■ days. It is well known for example (cf. and ref. therein) that if the 
^ ; coefficients of SDE are Liphitzian then the SDE generates a flow of homeo- 
morphisms, if coefficients are of the class C^^^ then SDE generates C"^-flow 
0^ '. of diffeomorphisms, equations for derivatives are obtained by formal differ- 

r ■ entiation of the SDE etc. 
^ ; Note that the similar questions for SDEs with reflection is much harder 

to answer. Even the problems about coalescence of two reflecting Brownian 
motions P, E], IH, [5] or differentiability of the Brownian reflecting flow (cr(x) = 
> ■ const) P, [7] need accurate and non-trivial considerations. 

The article below was published in Reports of Ukrainian Nat. Acad, of Sci. 
[B] (2005). Only some new references or minor remarks are added. 



(3 ! Assume that functions ak '■ satisfy the Lipschitz condition. Here 

■ = R^^^ X [0, oo). Consider an SDE in with normal reflection from the 

■ boundary: 

_> : 

"X ■ {dLpt{x) = ao{(pt{x))dt + YJk^i ak{^t{x))dwk{t) + 



{ +n^(dt,x), t G [0,r], (1) 

?o(x) = x, <^(0,x) = 0, X G R^, 

where {wk{t)^ A; = 1, . . . , m} are independent Wiener processes, n = (0, . . . , 0, 1) 
is a normal to hyperplane W^^^ x {0}, for each flxed x G a process ^(t, x) 
is non-decreasing in and 

^{t,x)= I I[{<^,(a;)eMd-ix{o}}'^(«^S) a:), 
JO 

i.e. ^(t, x) is increasing only on those instants of time when ift{x) G W^-^ X 
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{0}. Lipschitz property of the coefficients ensures ttie existence and tfie 
uniqueness to tfie solution of (1), cf. [S]. 

1. Existence of continuous modification. 

Theorem 1 [llOj]. There exists a modification of the processes (pt{x)^^{t^ x) 
(it will he denoted in the same way) such that 

1) for any x G M^l, the pair (ipt{x),^{t,x)), t >0, is a solution of (1); 

2) for any cj G H processes (pt{x), ^{t, x) are continuous in a pair of argu- 
ments (t, x), t > 0, x G M^. 

Tfie Tlieorem 1 is proved in a way similar to the corresponding proof 
used for the solution of SDE without reflection, cf. with the use of 
Kolmogorov's theorem on existence of continuous modification. 

It will be assumed further that cpt{x), ^{t, x) are already continuous. 

2. The joint motion of solutions started from different initial 
points. It is well known P that a solution of an SDE (without reflection) 
generates a flow of diffeomorphisms. However, the injectivity for reflecting 
flow can be failed as the following example shows. 

Example 1. Let d = 1, m = 1, ao = 0, ai = 1, i.e. (pt{x) is the reflected 
Brownian motion in Mi started from a: > : 



where t{x) is a moment, when the process x + w{t) gets zero for the first 
time. 

In other words, <-pt{x) is moving as x + w{t) before hitting 0, and then a 
motion of (pt{x) coincides with the reflected Brownian motion (pt{0) started 
from zero. 

The similar situation takes place in multi-dimensional space. 



<-Pt{x) = X + w{t) + ^(t, x), X > 0. 
It is easy to see that (pt{x)^^{t^ x) is of the form 
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Theorem 2. [TO] Denote by r{x) = M {t > : ipt{x) G R^"^ x {0}} the 

moment of the first hitting the hyperplane W^^^ x {0} by a solution started 
from X G R^. 

Then there exists a set Qq of probability 1 such that for all cj G the 
following statements hold true: 

1 ) for all x,y G R^, x y and t < mdux{r{x),r{y)} the inequality (pt{x) 
<-Pt{y) is satisfied; 

2) for any x G M.j_ there exists y = y{x,u) G R^^-^ x {0}, such that 
^t{x){x) = ^T{x){y) ifr{x) < oo. Moreover, 

(ptix) = (ptiy) fort > t{x). 

Remark. Informally this theorem can be formulated in the following way. 
A particle started from a point x G R^^^ x (0, cxd) does not hit any other 
particle before getting the hyperplane R*^^^ x {0}. At the instant t{x) it 
coalesces with some other particle, which started from R'^"^ x {0}. After this 
both particles moves together. 

3. Characterization of inner and boundary points of random set 

Theorem 3.fll| For almost all cu and all t G [0, T] the following equality 
of random sets takes place 

where t{x) = inf{s > : ips{x) G R^-i X {0}} is the moment of the first 
hitting the hyperplane R^^^ x {0} by the solution started from x. 

Moreover, for all R > Hausdorff measure H'^^ of the setdipt{R%)n{x G 
R% : \\x\\ < R) IS finite. 

4. Differentiability with respect to initial condition. 

As in Example 1, there is no reasons to expect that a solution of (1) is 
continuously differentiable in x even when coefficients of the SDE are infinite 
differentiable. However, it can be proved that for any t the mapping x 

(pt{x) belongs to a Sobolev space fl (R^,R'^). 

p>i ^' 

The equations for V(pt are not classical equations of stochastic analysis. 
We need the following definition. 
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Definition 1. Let wi{t), . . . ,Wm{t) be independent Wiener processes, 
jr^ = a{wkis), k = T~m, s < t), : x ^ M^, 6^ : x M^^ ^ W, k = 
0, . . . , m, and Xt be a continuous JT^-measurable stochastic process. Consider 
a random measure- valued process ^{t) = (5oI[{a;(i)=o}, where Sq is a probability 
measure on E, assigned unit mass to a point zero. 

A pair [yt, Zt) of jT^-adapted processes satisfies the equation 



dyt = ao{yt, Zt)dt + YJk=i ^k{yu Zt)dwk{t) - yt-du{t), 
dzt = bo{yt, zt)dt + YJk=i hivu zt)dwk{t),t > 0, 



(2) 



if: 

1) yt,t >0 has cadlag trajectories; 

2) Zt,t > has continuous trajectories; 

3) zt = zo + Jq bo{ys, Zs)ds + Y.T=i lo ^kiVs: Zs)dwkis),t > a.s.; 

4) for almost all uj the set {t > : Xt = 0} is contained in {t > : yt = 0}; 

5) for any stopping time r such that Xr ^ a.s., the following equality 
holds true 

ao{ys,Zs)ds + ^ / ak{ys, Zs)dwk{s) 
for all t e [t.t), T = in[{t > r : Xt = 0} . 



Theorem 4. fl2] Assume that functions ak-, hk-, k = 0,m satisfy Lipschitz 
condition. Then there exists a unique solution of (2) for any non-random 
initial condition (?/o,^o)- 

Theorem 5. /. If functions ak W^,k = 0, m satisfy 

Lipschitz condition then for a. a. uj a mapping 3 u ^ ^t{u) £ R"^ belongs 

to the space D W\^m^,,R'^) for a.a. t > 0. 

p>i ^' 

//. [12] Assume that functions ak '■ R+ W^,k = 0, m, are contin- 
uously differentiable and their derivatives are bounded. Suppose also that 
X^fcLi(<^fc,rf(^))^ > for all X G R'^^-^ x {0}, where ak^d is the d-th coordinate 
of a function ak = (afc,i, . . . , ak,d)^. 
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Then the Sobolev derivative VLpt{x) satisfies the SDE 



f m 

dVipt{x) = Vao{ipt{x))\/ipt{x)dt+ J2 ^ak{(pt{x))Vipt{x)dwk{t)- 

-PVipt-{x)n{dt,x), 
V(po{x) = 2, 

where 2 is an identity matrix, P is a matrix corresponding to the orthoprojec- 
tion on the d-th coordinate of the space M.^, n{dt, x) is a point random measure 
such that n{{t},x) = 1 iff ipt{x) belongs to the hyperplane W^^^ x {0}. 

Remark. Equation (3) is understood in the sense of Definition 1. In this 
case we take the d-th coordinate of <^t{x) as x^, the d-th row of V(pt{x) as yt, 
the first {d — 1) rows of Vcpt{x) and <-pt{x) as the process Zt. 

Remark. The process V(ft{x) can be chosen measurable in t^x^co. 

Remark. The similar result for constant diffusion coefficient was obtained 
in |[7j. Moreover, it was proved that for all x and a. a. to the mapping ipt is 
continuously differentiable in some neighborhood of x. 

Let us compare Sobolev differentiability and usual differentiability of the 
mapping (ft{-,cj). It is well known that if the diffusion matrix is a constant 
then for a.a. u and all t the mapping x Lpt{x^Lo) satisfies Lipschitz con- 
dition. Therefore (pt{x) is differentiable for A -a.a. X G G by Rademacher's 
theorem [H]. Since the usual and Sobolev derivatives are equal (if they ex- 



ist), so equation for usual derivative coincides with that for Sobolev. It is 
not difficult to prove that the usual derivatives exist not only for 
cj, but for all x and a.a. to. However almost everywhere local continuous 
differentiability is not evident. 

It should be noted that [7] does not imply that for a.a. uj the mapping 
X (pt{x,co) is continuously differentiable. Really, for the process from 
Example 1: 



P{x ^ <-Pt{x) is continuously differentiable) = 

but for each xq > : 
P{x ^ <-Pt{x) is continuously differentiable in some neighborhood of xq) = 1! 

The fact that 3 x ^ ^t{x) is not continuously differentiable seems to 
be typical, because lankV (pt{x) < d — 1 if r{x) < t and i-Pt{x), r{x) > 
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coincides with the flow without reflection, so det V(/?i(a:), r{x) > \\x\\ < r, 
is separated from zero for any r > 0. 

5. Absolute continuity of image-measures driven by (pt. Let be a 
finite measure in which is absolute continuous w.r.t. Lebesgue measure. 
Consider a measure- valued process = /i o (^^\ ^ > 0. 

Let us introduce a random set Ot{io) = {x G M.j_ : t < r{x)}, where 
t{x) = inf{s > : (ps{x) G R^~^ x {0}} is the moment of the first hitting the 
hyperplane W^~^ x {0} by the process (ps{x). 

Theorem 6. |13j For a. a. to and every t > a measure fit is represented 
as a sum of orthogonal measures fit = /ij^ o (^^-^ + fJi\^d\^Q o ^t^-, such that 

a) the first measure is absolute continuous w.r.t. d-dimensional Lebesgue 
measure and the second one is singular; 

b) the support of measure (j\^d\^Q °^t^^ is contained in the set (pt{M.^^^ x {0}) 

of the a-finite {d — 1) -dimensional Hausdorff measure H'^~^. 

The proof of the first part of the theorem follows from [TB] , and the second 
part follows from Theorem 3. 

The next theorem gives a sufficient condition that ensures the absolute 
continuity of ^^t'^ with respect to i/^'^j^^ ^^^^y which is the restriction 

of F^-i to the set dLpt{M.%). 

Theorem 7. fl3] Assume that for fi-a.a. x G : 

P(rank V(/?f(x) >d-l,t>Q) = l. (4) 

Then with probability 1 the absolute continuity 

holds for all t > 0. 

Here f^\^d\^Q , ^^^^Ig^ (K<i) ^^^^ restrictions of measures ^u, H^~^ to the sets 

W^\Ot, d(pt{Wl_), respectively. 

Remark. In contrast to Hausdorff measure H'^^^ in M.^, its restriction to 
the set diptiWl^) = (ft(^'i \ Ot) is a cr-finite measure (Theorem 3). So the 
notion of absolute continuity in (5) does not require any specification. 

The proof is provided by using the co-area formula [H] similarly to the 
case m = n, cf. ^E\, [151 
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The verification of Theorem 7 conditions is quite difficult. Let us give 
more simple sufficient conditions that ensure (4). Assume that functions 
ak,0 < k < m, are continuously differentiable. 

Denote by Ust{x), s <t, a solution of the following linear SDE: 

dUst{x) = VaQ[Lpt{x))Ust{x)dt + YJk=i^^k{^t{x))Ust{x)dwk{t),t > s, 

Let us represent a random set A = At{x) = {s G [0,t] : (pi{x) > 0} as 
a disjoint union of random intervals A = [ao{x) , Pq{x)) U {ai{x) , Pi{x)] U 

oo 

IJ {ak{x) , Pk{x)) , where ao{x) = 0,Pi{x) = t. Let P be the same as in The- 

k=2 

orem 5. 

Theorem 8. fl3] Let the conditions of the second part of Theorem 5 be 
satisfied, t > 0. Assume that for all k > and a. a. x G U 

P (rank((l - P)U^,p,{x){l - P)) = d - 1) = 1. (6) 

Then relation (4) holds true. 

Corollary. Assume that for all x G W^, s < t : 

P (^(ffix) = 0, det =0^ =0, where Ust{x) is the matrix getting out 

from Ust{x) by deleting last raw and last column. Then (6) fulfills. 

Observe that the assumption of the Corollary is the requirement of non- 
hitting zero by two-dimensional Ito process, and this is usually easier to check 
than (6). 

Example 2. Let d = 2, then assumptions of the Corollary to the Theorem 
8 are satisfied if for all 7^ vectors {ak^i{x))i<k<m and 

{V yak^2{x))i<k<m, are linear independent. 

6. Flows generated by SDE with reflection in arbitrary set. Let 

G C be a closed set with smooth boundary such that the following SDE 
with normal refiection on the boundary of G has a unique strong solution 
defined for all x G G, i > : 

d(pt{x) = ao{(pt{x))dt + XlfcLi ak{(pt{x))dwk{t)^ 

-\-n{ipt{x))^{dt,x), t>0, X gG, 
cpo{x) = X, ^{0, x) = 0, ^{t, x) = /q 'S.{^,{x)edG}^{ds, x), 
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where n{x) is the inward normal at a boundary point x G dG, ^{t^x) is 
continuous and non-decreasing in t process for every fixed x E G. 

Assume that for any x G dG there exists a neighborhood 0{x) and G^- 
diffeomorphism which transforms the set 0{x) fl G into {x G : < 
1, > 0} in such a manner that Vax{y) = n{y) = (0, . . . , 0, 1), y G dG fl 
0{x). 

Applying a localization of solutions, all statements on fiows in half-space 
can be easily generalized (of course with natural changes) to the case of SDE 
in G. For example, relation (6) will be of the form: 

P (rank(P(^^Jt/«,^,(x)P(^«J) = d - 1) = 1, (7) 

where P{x) is orthoprojection on the orthogonal complement to ri{x) where 
X G dG. 

Example 3. Let G = {x G : < 1} be a unit disk, ipt{x),t >0,xE 
G be a Brownian motion in G with refiection on the boundary. I.e. (pt{x) is 
a solution of the SDE 

dipt{x) = dw{t) + n{ipt{x))^{dt, x), t > 0, 
(po{x) = X, ^(0, x) = 0, X G G, 

where w{t) is a two-dimensional Wiener process. 

Let us describe inner and boundary points of the set (pt{G). Now, the 
stopping time r(x) from Theorem 2 is of the form r(x) = inf{t > : 
X + w{t) G dG}. So, the set of inner points of '■Pt{G) is equal to {x + w{t) : 
X E G^t < r(x)}. 

Introduce a stopping time a = mf{t > : ||w^(^)|| = 2}. Observe that 
sup^; t(x) < a. Therefore the analogues of Theorems 2,3 imply that for every 
t > a the random set <^t{G) coincides with the nowhere dense set (pt{dG) of 
finite Hausdorff measure H^. Moreover, 

P (Vt > (7 Vx, ||x|| <13yedG,xj^y: iptix) = ^t{y)) = 1- 

It is interesting to compare this result with \T\, where it is proved that any 
two solutions of (8) started from different initial points of G never meet each 
other with probability 1, that is 

yx,yeG,x^y: P (3t > : ipt{x) = ipt{y)) = 0. 
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Note that in this example Ust{x) is an identity matrix, so condition (7) 
is obviously satisfied. Thus for any absolute continuous measure /i on G we 
have the absolute continuity 

with probability 1. In particular, for a. a. u and all t > a: 

Observe that if G is not a unit disk but any domain with "nice" boundary, 
and the boundary does not contain any two perpendicular segments, then (7) 
and so (4), (9) are also satisfied. Exactly the same condition on the boundary 
appears in [H]. Moreover, it can be easily shown that if the boundary contains 
two perpendicular segments then (4), (7), and (9) are false. 
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